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Abstract 



We give a formulation of the single particle occupation probabilities for a 
system of identical particles obeying fractional exclusion statistics of Haldane. 
We first derive a set of constraints using an exactly solvable model which 
describes an ideal exclusion statistics system and deduce the general counting 
rules for occupancy of states obeyed by these particles. We show that the 
problem of negative probabilities may be avoided with these new counting 
rules. 
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I. INTRODUCTION 



A few years ago Haldane proposed a generalized exclusion principle in which adding 
particles in to a system leads to a change in the dimension of the single particle space. 
Specifically, the generalized exclusion principle envisages systems in which the addition of 
one particle blocks off g single particle states for the others, where g, the fractional exclusion 
statistics(FES) parameter is an arbitary number. Obviously g = for bosons and g = 1 for 
fermions. This leads to the following formula for -Djv(g, d), the dimension of the N particle 
Hilbert space, if the dimension of the single particle space is d, 

(d+(l-g)(N-l)\ 
DnM = N\(d-l-g(N-l)\ (1) 

which reduces to the familiar expressions for Bose and Fermi statistics for g = and g — 1. 

Thermodynamic properties of an ideal gas of exclusion particles have been investigated 
widely. Specifically, a definition of an ideal gas of particles with nontrivial exclusion statistics 
was given in references 0JJ]. In this definition it was assumed that if there were d levels of 
energy e, then the dimension of the Hilbert space with N particles of energy e is given by 
equation (fj). The distribution function may then be computed and is given by 

<e) = 1 , (2) 
(w(e)+g) 

where w(e) is the solution of the equation 

w(e) 9 (l + w(e)) {1 ~ 9) = eP* (3) 

and f3 is as usual the inverse temperature. If we attempt to interpret this distribution 
function as arising from the statistical mechanics of a single mode with the statistical weight 
p n &~ for the mode to be occupied by n particles, then it was found |||7] that some of 
the p n s are invariably negative if g is different from and 1. This raises the problem of 
interpreting these negative probabilities. It has also been speculated that these negative 
probabilities are an essential feature of nontrival exclusion statistics M. 
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A crucial property of exclusion statistical interactions is that they should cause shifts 
in single particle energies at all scales ||. This property is realized by a large class of one 
dimensional models of interacting fermions where Fermi liquid theory breaks down PJT0|] . In 



fact it has been shown exactly that quasiparticles with nontrival exclusion statistics exist in a 
class of models that are solved by the Bethe ansatz [0,|nj] . In particular the quasiparticles of 



the Calogero- Sutherland model (CSM) behave like ideal exclusion statistics system PJTTHIBl. 
A feature of the exclusion statistics as gleaned from the analysis of various models is that 
the exclusion acts across a set of levels unlike in the case of Fermi or Bose statistics where 
the exclusion principle is stated with a single level in mind. It is this crucial difference 
that results in the occurance of negative probabilities. We will show that the particles 
obeying fractional exclusion statistics can be characterised by constraints on the sets of 
occupation numbers. There are no negative probabilities if these constraints are obeyed. 
If these constraints are relaxed then the negative weights arise in order to compensate for 
the resulting over counting. Indeed this is the way we encounter negative probabilities in 
other systems in physics- for example in gauge theories, they arise in the ghost sectors. 
Ghosts come from the Jacobian associated with nonlinear gauges which essentially ensure 
the correct counting of states. Another example is that of the Wigner distribution function 
in quantum mechanics which is not positive definite precisely because some constraints are 
relaxed. A formulation based on the variable number of single particle states, which depends 



on the total number of particles in the system, has been discussed by Isakov [14]] as a way to 



avoid the problem of negative weights. Recently, a microscopic interpretation of exclusion 



statistics systems has been advanced by Chaturvedi and Srinivasan [|15|] where they show 
how this problem of negative probabilities may be solved for semions, g = 1/2. They have 
also indicated how their method may be generalised to other values of g. 

In this paper we first discuss the origin of negative probabilities in exclusion statistics 
particle systems. To do this we have chosen an unusual starting point in an equation and its 
solution given by Ramanujan (jnj. This starting point makes precise the statements about 
the occurence of negative probabilities. We then formulate a counting principle based on 



the set of constraints which reproduces the Haldane dimension formula. We first extract 
the counting rules starting from an exactly solvable model of interacting particles and state 
them in the form of counting rules for arbitrary systems obeying exclusion statistics. This 
method not only avoids the negative probabilities, but with minimal modification reproduces 
the results derived by Chaturvedi and Srinivasan for the semion. The counting principle is 
however not restricted to semions alone. 



II. THE PROBLEM OF NEGATIVE WEIGHTS 

The problem of negative probabilities was first pointed out by Nayak and Wilczek || and 
elaborated by Polychronokos in a recent paper 0. In order to clarify the origin of negative 



probabilities or weights, we first discuss an equation and its solution due to Ramanujan p6 
Ramanujan considered the following equation: 

aqX p - X q + 1 = 0, (4) 

where a may be complex and p, q are positive. The general solution for X d is, 

oo 

X d = J2C N (p,q,d)a N , (5) 

N=0 

where C (p, q, d) — 1 and Ci(p, q,d) = d and 

d N ~ x 

C N (p,q,d) = — XKd + Np-jq), N>2. 

To make connection with the result obtained by Polychronokos [|7] , which is a particular 
case of the general solution given by Ramanujan, we now put p — (1 — g) and q — 1, then 

C ^~9,hd) = d m{d _ gN)] (6) 

which is clearly different from the dimension formula of Haldane. However, it correctly 
reproduces the bosonic and fermionic dimension formula for g = and g = 1 respectively. 
This dimension formula was derived independently by Polychronokos with the restriction 



that any two particles are atleast g sites apart when placed on a periodic lattice. One can 
also derive the Haldane dimension formula with the restriction that any two particles are g 
sites apart but without the restriction of periodicity. 

Further if we put X = (1 + w; -1 ) and a = e~^ e in eq.(^), we immediately obtain equation 
(0) derived earlier by Wu. The important point to notice here is that the dimension formula 
that precisely leads to the distribution function derived earlier |2|-|5| is given by Cat and not 
the Haldane dimension formula. In the limit d » 1, however, it is easy to see that 

C N (l-g,l,d) = D N (g,d) + 0(^). 

Therefore in the continuum limit, the Cn and Dn are approximately the same. 
The grand canonical partition function of the system may be written as, 

oo 

Z = (l + w-'f = Y,C N (l-g,l,d)e-P N % (7) 

N=0 

where w satisfies equation (|3|). We have also assumed that all the energy levels are degenerate 
with energy given by e. Note that this is an exact expression and no assumption is required 
on the single particle dimension d. The negative weights arise [§,0, when one insists on 
expanding 1 + w' 1 in powers of e _/3< \ From equation (|5J) and the definitions following the 
equation, it follows that, 

oo 

1+uT 1 = £C B (l-0,l,l)e-* M (8) 



n=0 



The weights 



0,(1-2,1,1)= Pn = na--) w 



m=2 771 



are always negative for gn > m for some m J/J. This is indeed the problem of negative 
weights associated with exclusion statistics and is claimed to be inherent in the exclusion 
statistics. There are however a few points to note: The negative probabilities arise because of 
our insistence on the factorization [15] implied in eq.([7|). For example, combining equations 



(0) and (H) we have, 



z = E (11^(1-^1,1))^^^, (io) 

K> 3 

where the sum is an unconstrained one over all sets of occupation numbers. The over 
counting resulting from this unconstrained sum is compensated by the occurence of negative 
weights. We next derive the precise counting rules which impose constraints on this sum 
and avoids this problem. 



III. REALIZATION IN AN INTERACTING SYSTEM AND COUNTING RULES 

Any realization of fractional exclusion statistics must have its origins in systems of in- 
teracting particles. The expectation is that under certain conditions systems of interacting 
particles which obey Fermi or Bose statistics may be described in terms of quasiparticles 
(or quasiholes) which obey fractional statistics. The quasiparticles of the CSM behave like 
ideal exclusion statistics particles. The main feature of CSM is that the total energy of the 
many-body system can be written in terms of single quasi-particle energies which involve 
shifted momenta and these shifts contain the information about the exclusion statistics of 
the quasiparticles. In this section we analyse these shifted momenta and make explicit con- 
nection with the formula in eq. ([!]). We then use them to obtain constraints on the allowed 
set of occupation numbers. These are what we refer to as the counting rules that reproduce 
the formula in eq.([TJ). The statistical mechanics of the system obeying these constraints is 
then the same as that defined by Wu || and all statistical weights are positive. 



We begin with the trignometric Sutherland model |17| of an N-particle system on a ring 
of unit radius. The Hamiltonian is given by, 

H _ f d 2 , y 2g(g-l) 

kdxrfysm 2 [( Xi - Xj )/2y [ } 

where g is the interaction parameter. We will soon identify this with the statistical parameter 
of the exclusion statistics. While the model can be applied to both interacting bosons and 
fermions, we choose to work in the fermionic basis here after. The energy of an N-fermion 
state may be written in terms of shifted momenta as 



E^tfm, (12) 



where = 0, 1 and the shifted momenta ki (also called pseudo momenta in Ref. Jy|) are 
given by 

h = m-il-g) ^* ~ N ^\ (13) 

where va\ are distinct integers, N^ + ' are the number of particles with shifted momenta less 
(greater) than ki. Note that we could have also started with the Calogero- Sutherland model 
with harmonic confinement. The results below follow analogously with the proviso that we 
have shifted energies instead of shifted momenta. 

First we establish the relationship between the shifted momenta given above and the 
Haldane's dimension formula (p]). Consider the above system with an upper and lower 
cutoff on the momenta, k max and k m i n respectively. We divide this range of momenta into 
cells of unit length (the first and last cells could be smaller) and define the occupancy of 
the j th cell, rij to be the number of particles with momenta ki, such that j + 1 > ki > j. 
We identify single particle space dimension d with the number of cells in the range, i.e 
d = k max — k min , where d may be fractional. If we now denote the range of the m^s by dp, 
we have 

d F = m max -m min = d+ (1 - g)(N - 1). (14) 

Since there exists an for every ki, the total number of states in the range k max — k m i n is 
the same as that between m max — m m i n . The total number of states is then the number of 
ways N distinct integers can be picked from d F distinct integers, i.e dp CN, as in fermionic 
description. Substituting for dp from the above expression we immediately reproduce the 
Haldane dimension formula in eq.([l]) ||. 

In order to obtain the counting rules we will first derive three properties of the set of 
momenta {ki}. If ki are ordered such that they increase with increasing i, then we have, 
k i+ i — ki = m i+ i — rrii — (1 — g). If g < 1, then it follows that m i+ i > Further, if 



rrii + i — rrii = 1 then ki + \ — ki = g and if raj+i — mi > 1 then fcj+i — ki > (1 + g) > 1 because 
mj's are integers. 

We can then draw the following three conclusions from the properties of the shifted 
momenta kf. 

1. The ordering in fc^s is the same as the ordering in m^s. 

2. "Close packed" m^s with unit spacing correspond to "close packed" kiS with spacing 
9- 

3. The gaps between any two non-close packed /c,s is greater than 1. Therefore all the 
kiS in any cell are close packed. 

We now come to the question as to what are the constraints on the sets of occupation 
numbers {ki}. For example, if g = 0, there are no constraints as in the bosonic case. If 
g = 1 the constraints are rij < 1 as in the case of fermions. For any other g, one obvious 
constraint come from the second property derived above, namely the occupancy of the j-th 
cell rij < i which specifies the maximum occupancy of a given cell assumed to be of unit 
spacing. This is the same constraint one derives from the distribution function of Wu (0). 
An important departure from the usual bosonic and fermionic case is that the cell size is 
important and cannot be arbitrarily taken to zero as in the case of bosons and fermions || . 

There are further constraints on the occupancy. To formulate them we use the third 
property. Let k^ be the lowest momentum in the j th cell. Then from the second and third 
property, it follows that 

k L + g{n j -l)<j + l. (15) 

We can write ki as ki = j + f{ki), where f{ki) denotes the fractional part of k^ that is, 
< f(k L ) < 1. We then have, 

f(k L )+g{ nj -l) < 1. (16) 

From equation (|13D , we can express /(^l) as a function of the occupation numbers, 
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f -(l-g) 



[N-j - jg - ( Wj - - 1)] 



(17) 
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where N~ = J2i<j n i an d = J2i>j n i- Equations fllED and ( |TT| ) then constitute a set of 
constraints on the occupation numbers. 

We will now show that these form a complete set of constraints. Namely, given any set 
of occupation numbers, {nj}, that satisfies the constraints, there exists a set of momenta, 
{h}, that realizes it. To do this, consider a set {%}, where j m i n < j < j ma x- The lowest 
value of the momentum in the j th cell is uniquely determined by the occupation numbers 
through equation([T7|). Because of the third property, all the other momenta are also uniquely 
determined. Hence we have shown that there are no more constraints. Equations fllCf ) and 
(|17|) form a complete set of constraints. Note also that the above logic implies that there is 
a one to one correspondence between the sets of occupation numbers, {nj}, that satisfy the 
constraints ( |i6| ) and (|H])and the sets of momenta, {h}, that satisfy equation (Jl3|). 

We can now remove the scaffolding of the Sutherland Model that we started with and 
define exclusion statistics system by the above constraints. The connection to the dimension 
formula in eq.(0) established earlier implies that 



where N = J2j n j an d -F({n,}) = 1 if {nj} satisfy the constraints and zero otherwise. Note 
that the weights now are positive definite. There are no negative weights once the constraints 
are imposed. 

Next, we construct some simple examples from the above counting rules. For simplicity 
we look at occupation numbers for special values of g = 1/m where m is an integer. The 
rules formulated above for the occupation number of exclusion particles may be combined 
and restated thus: 

Let m = 1/g, and let Ni be the number of particles in the occupied states below 
some ith level, Ni = J2j<i n i- Then an occupation ni(rii < m) is allowed iff' 
(Ni mod m) < (m — n«) . 



= D N (g,d) 



(18) 



{ni} 
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This rule now includes all the three constraints stated above. 

To see how this rule is implemented, consider a system of N-particles spread over d 
states. In order that D n is an integer, we choose N = mp + 1, where p is an integer. Since 
iV < md, we have p < d. We shall divide these d states into cells. An allowed configuration 
may be represented as a string of numbers (n^r^, ...), where each rii < m denotes the 
occupancy of levels ordered from left to right. Instead of dealing with a configuration where 
all N particles are spread over d states (some which may be empty), we can simplify the 
discussion by considering one cell at a time. Each cell may now have a partition of m. This 
allows us to fill the subsequent cells without reference to the previous cell according to the 
counting rules since JVj mod m — 0. We now fill each cell with a partition of m which is 
allowed by the rules given above. This then generates all possible allowed configurations 
whose sum is given by D N . 

If, in particular, we are interested in expectation values of symmetric functions of rii, 
we can work with symmetrised weights. Consider a symmetric operator 0({nj}). The 
expectation value of this operator may be written as, 

<0(W)>=?f™, (19) 

where 

F s ({ni}) = ^X>W. (20) 

Here p stands for all permutations of the allowed configurations. Every allowed configuration 
in {rii} may be characterised by the multiplicities q n , namely a given allowed configuration 
may be written as a string, m qm (m — l) qm - 1 ...l qi , where q 1 + 2q 2 + ... + mq m = N. We 
may now also allow any permutation of these occupancies (with zeros added to make up 
d-states). The dimension of the N-particle space may then be written as 

D N (g,d)=J2f»(q u q 2 ,...) d C q , (21) 

{qn} 

where q = J2™=i In- The new weights / are defined as, 
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eN / \ M a (q 1 ,q 2 , ...q m ) / 00 s 
f m 92, -q m ) = TTT— ~ y ( 22 ) 

where M a allowed configurations after symmetrising and M t is the total number of config- 
urations for a given set of g's which define a configuration. We shall clarify this now with 
specific examples. 

A. The case of semion (g = 1/2, m = 2) 

The maximal occupancy of a state in this case is 2. Hence allowed occupancy of a state 
is 2 or 1. Zeros may occur any where without changing the rules. Let us implement this in 
the specific case of d = 4, N = 5, say. In this case the allowed configurations are given by the 
strings (2210), (2111), (1121). In the first configuration, zero can be anywhere and therefore 
there are four configurations. Notice that a string of the form (1211) or (1112) violates the 
counting rules. Therefore counting all the allowed configurations we obtain D 5 (l/2,4) = 6. 
This is exactly what one gets from the Haldane formula. 

Further if we symmetrise each of these allowed configurations, then the new weights may 
be computed using eq.(^). In the specific case of m = 2, we have 

M t {qi, <&) = qi+q2 C q2 , M a ( qi , q 2 ) = p C q2 , (23) 

where p is defined through the equation N = 2p+l. The corresponding / is therefore given 
by 

/" = (24) 

Note that these weights, wheather in the symmetrised form or unsymmetrised form, are 
positive definite. Further, this is exactly the formula derived by Chaturvedi and Srinivasan 
|T5| in their microscopic analysis of Haldane statistics for semions. 

It is important to stress the differences in these two approaches- in their analysis 
Chaturvedi and Srinivasan start from a formulation of the statistical mechanics of a system 
by removing factorizability of the weights as a criterion. They derive the expression for 
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the weights in eq.(|24|) by imposing the conditions positivity and the requirement of sym- 
metry (all configuratons which are permutations of each other carry the same weight). Our 
starting point is the Sutherland model. We derive our rules from the properties of shifted 
momenta. After removing this scaffolding, we obtain not only positive definite weights for 
each configuration but when symmetrised they reproduce the results of Chaturvedi and 
Srinivasan. 



B. The case with g = 1/3 or m = 3 

The maximal occupancy of a state in this case is 3. The allowed configurations for each 
cell are (3), (21), (12), (111). That is we can form a string of allowed configuration with any 
of these cells in any order to make up N particles. Any number of zeros may be added in 
between to make up a total of ci-states. 

As in the semion case we may consider expectation values of symmetric functions of rij. 
Following the same procedure we can derive the symmetrized weights defined in eq.(p^). 
Since, m = 3, we have 

M t ( qi ,q 2j q 3 )=^ + ^C q3 ^C q2} M a (q u q 2} q 3 )=^C q3 ™C q2 (2f\ (25) 

where p, as before, is defined through the equation N = 3p + 1. The corresponding weight 
/ is therefore given by, 

f N = 93 Q2 ± ' (26) 

, ' 3 gi+<22+<?3f7 <?i+<32f7 ^ ' 

°93 ^92 

These weights are again positive definite. Chaturvedi and Srinivasan also suggest how 
their method may be extended beyond the semion case which they considered in detail. 
However, this extension requires additional conditions which are not imposed in the semion 
case. In contrast, our rules as derived from the point of view of an exactly solvable model 
are completely specified independent of the actual value of g (or m) . There is an algorithm 
to derive for arbitrary m though this gets complicated for larger m. 
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IV. SUMMARY 



To summarise, we have analyzed the origin of negative probabilities in exclusion statis- 
tics systems. To do this we have chosen an unusual starting point in an equation and its 
solution given by Ramanujan. This starting point makes precise the statements about the 
occurence of negative probabilities. Further, we have formulated a counting principle which 
reproduces the Haldane dimension formula. It can therefore be used to define exclusion 
statistics purely in terms of state counting. The negative probabilities discussed in liter- 
ature can be understood as arising when the system constrained by the counting rules is 
replaced by an unconstrained one. The negative weights then compensate for the introduc- 
tion of unphysical configurations. This is therefore exactly analogous to other situations in 
physics where negative probablities arise, for example, the ghosts and negative norm states 
in gauge theories or as in the case of Wigner distribution in quantum mechanics. 
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